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Abstract 

The dynamics of the equal-time cross-correlation matrix of multivariate financial 
time series is explored by examination of the eigenvalue spectrum over sliding time 
windows. Empirical results for the S&P 500 and the Dow Jones Euro Stoxx 50 
indices reveal that the dynamics of the small eigenvalues of the cross-correlation 
matrix, over these time windows, oppose those of the largest eigenvalue. This be- 
haviour is shown to be independent of the size of the time window and the number 
of stocks examined. 

A basic one-factor model is then proposed, which captures the main dynamical 
features of the eigenvalue spectrum of the empirical data. Through the addition of 
perturbations to the one-factor model, (leading to a 'market plus sectors' model), 
additional sectoral features are added, resulting in an Inverse Participation Ra- 
tio comparable to that found for empirical data. By partitioning the eigenvalue 
time series, we then show that negative index returns, (drawdowns), are associated 
with periods where the largest eigenvalue is greatest, while positive index returns, 
(drawups), are associated with periods where the largest eigenvalue is smallest. The 
study of correlation dynamics provides some insight on the collective behaviour of 
traders with varying strategies. 
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1 Introduction 



In recent years, the analysis of the equal-time cross-correlation matrix for a va- 
riety of multivariate data sets such as financial data, [Tf2 3 4 5 6 7 8 i9]fl0llll|ll2] .[ 
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electroencephalographic (EEG) recordings, [TSUI], magnetoencephalographic 
(MEG) recordings, [15], and others, has been studied extensively. Other au- 
thors have investigated the relationship between stock price changes and liq- 
uidity or trading volume, [H)jl7|ll8] . In particular, Random Matrix Theory, 
(RMT), has been applied to filter the relevant information from the statistical 
fluctuations, inherent in empirical cross-correlation matrices, for various finan- 
cial time series, [Tf2 ,3 4 5 6 7 8 9|TT0] . By comparing the eigenvalue spectrum of 
the correlation matrix to the analytical results, obtained for random matrix 
ensembles, significant deviations from RMT eigenvalue predictions provide 
genuine information about the correlation structure of the system. This infor- 
mation has been used to reduce the difference between predicted and realised 
risk of different portfolios. 

Several authors have suggested recently that some real correlation information 
may be hidden in the RMT defined random part of the eigenvalue spectrum. 
A technique, involving the use of power mapping to identify and estimate the 
noise in financial correlation matrices, [19], allows the suppression of those 
eigenvalues, associated with the noise, in order to reveal different correlation 
structures buried underneath. The relationship, between the eigenvalue den- 
sity c of the true correlation matrix, and that of the empirical correlation ma- 
trix C, was derived to show that correlations can be measured in the random 
part of the spectrum, [20|l2Tj . A Kolmogorov test was applied to demonstrate 
that the bulk of the spectrum is not in the Wishart RMT class, [22] , while the 
existence of factors, such as an overall market effect, firm size and industry 
type, is due to collective influence of the assets. More evidence that the RMT 
fit is not perfect was provided, [23], where it was shown that the dispersion 
of "noise" eigenvalues is inflated, indicating that the bulk of the eigenvalue 
spectrum contains correlations masked by measurement noise. 

The behaviour of the largest eigenvalue of a cross-correlation matrix for small 
windows of time, has been studied, [24J, for the DAX and Dow Jones Industrial 
average Indices (DJIA). Evidence of a time- dependence between 'drawdowns' 
('draw-ups') and an increase (decrease) in the largest eigenvalue was obtained, 
resulting in an increase of the information entrop^pi the system. Similar 
techniques were used, [25], to investigate the dynamics between the stocks of 
two different markets (DAX and DJIA). In this case, two distinct eigenvalues of 
the cross-correlation matrix emerged, corresponding to each of the markets. 
By adjusting for time-zone delays, the two eigenvalues were then shown to 
coincide, implying that one market leads the dynamics in the other. 

Equal-time cross-correlation matrices have been used, |26j, to characterise dy- 
namical changes in nonstationary multivariate time-series. It was specifically 



In information theory, the Shannon entropy or information entropy is a measure 
of the uncertainty associated with a random variable. 
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noted that, for increased synchronisation of k series within an M— dimensional 
multivariate time series, a repulsion between eigenstates of the correlation ma- 
trix results, in which k levels participate. Through the use of artifically-created 
time series with pre-defined correlation dynamics, it was demonstrated that 
there exist situations, where the relative change in eigenvalues from the lower 
edge of the spectrum is greater than that for the large eigenvalues, implying 
that information drawn from the smaller eigenvalues is highly relevant. 

The technique was subsequently applied to the eigenvalue spectrum of the 
equal time cross-correlation matrix of multivariate Epileptic Seizure time series 
and information on the correlation dynamics was found to be visible in both 
the lower and upper eigenstates. Further studies of the equal-time correlation 
matrix of EEG sig nals, |13j . investigated temporal dynamics of focal onset 
epileptic seizures^ and showed that zero-lag correlations between multichannel 
EEG signals tend to decrease during the first half of a seizure and increase 
gradually before the seizure ends. A further extension (to the case of Status 
Epilepticus, [14]), used the equal-time correlation matrix to assess neuronal 
synchronisation prior to seizure termination. 

Examples have demonstrated, [27], that information about cross correlations 
can be found in the RMT bulk of eigenvalues and that the information ex- 
tracted at the lower edge is statistically more significant than that extracted 
from the larger eigenvalues. The authors introduced a method of unfolding 
the eigenvalue level density, through the normalisation of each of the level dis- 
tances by its ensemble average, and used this to calculate the corresponding in- 
dividual nearest-neighbour distance. Those parts of the spectrum, dominated 
by noise, could be distinguished from those containing information about cor- 
relations. Application of this technique to multichannel EEG data showed the 
smallest eigenvalues to be more sensitive to detection of subtle changes in the 
brain dynamics than the largest. 

In this paper, we examine eigenvalue dynamics of the cross-correlation matrix 
for multivariate financial data with a view to characterising market behaviour. 
Methods are reviewed in Section [2j the data described in Section [3] and in 
Section H] we look at the results obtained both for the empirical correlation 
matrix and the model correlation matrices described. 



2 A partial or focal onset seizure occurs when the discharge starts in one area of 
the brain and then spreads over other areas. 
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2 Methods 



2. 1 Empirical Dynamics 



The equal-time cross-correlation matrix, between time series of equity returns, 
is calculated using a sliding window where the number of assets, N, is smaller 
than the window size T. Given returns Gi(t), i = 1,. . . , iV, of a collection 
of equities, we define a normalised return, within each window, in order to 
standardise the different equity volatilities. We normalise Gi with respect to 
its standard deviation <jj as follows: 

9l (t) = Gl (0 = Gi (0 (1) 



Where <7j is the standard deviation of Gi for assets i — 1, . . . , N and Gi is the 
time average of Gi over a time window of size T. 

Then, the equal-time cross-correlation matrix is expressed in terms of (t) 
C y = (f) ft- (f)> (2) 



The elements of Cij are limited to the domain —1 < < 1, where CV,- = 1 
defines perfect positive correlation, = —1 corresponds to perfect negative 
correlation and = corresponds to no correlation. In matrix notation, the 
correlation matrix can be expressed as 

C = ^GG T (3) 



Where G is an iV x T matrix with elements g it . 

The eigenvalues Aj and eigenvectors Vj of the correlation matrix C are found 
from the following 

Cv, = A,v, (4) 



The eigenvalues are then ordered according to size, such that Ai < A 2 < . . . < 
Aat. The sum of the diagonal elements of a matrix, (the Trace), must always 
remain constant under linear transformation. Thus, the sum of the eigenvalues 
must always equal the Trace of the original correlation matrix. Hence, if some 
eigenvalues increase then others must decrease, to compensate, and vice versa 
(Eigenvalue Repulsion). 
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There are two limiting cases for the distribution of the eigenvalues [T3f26] . 
When all of the time series are perfectly correlated, Ci ml, the largest eigen- 
value is maximised with a value equal to N, while for time series consisting of 
random numbers with average correlation Cj ~ 0, the corresponding eigenval- 
ues are distributed around 1, (where any deviation is due to spurious random 
correlations). 

For cases between these two extremes, the eigenvalues at the lower end of the 
spectrum can be much smaller than \ max . To study the dynamics of each of 
the eigenvalues using a sliding window, we normalise each eigenvalue in time 
using 

fat) - 

x ^ = V J ( 5 ) 



where A and a x are the mean and standard deviation of the eigenvalue % 
over a particular reference period, r. This normalisation allows us to visually 
compare eigenvalues at both ends of the spectrum, even if their magnitudes 
are significantly different. The reference period, used to calculate mean and 
standard deviation of the eigenvalue spectrum, can be chosen to be a low 
volatility sub-period, (which helps to enhance the visibility of high volatility 
periods), or the full time period studied. 



2.2 One-factor Model 



In the one-factor model of stock returns, we assume a global correlation with 
the cross-correlation between all stocks the same, po- The spectrum of the 
associated correlation matrix consists of only two values, a large eigenvalue 
of order (N — l)po + 1, associated with the market, and an (N — 1)— fold 
degenerate eigenvalue of size 1 — po < 1- Any deviation from these values is 
due to the finite length of time series used to calculate the correlations. In 
the limit N — > oo (even for small correlation, i.e. p — > 0) a large eigenvalue 
appears, which is associated with the eigenvector vi = (77^) (l,l,l...l), and 
which dominates the correlation structure of the system. 



2.3 'Market plus sectors' model 



To expand the above to a 'market plus sectors' model, we perturb a number of 
pairs iV of the correlations po + p n , where — 1 — po < Pn < 1 — Po- Additionally, 
we impose a constraint V] Pn = 0, ensuring that the average correlation of the 

N 
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system remains equal to po- These perturbations allow us to introduce groups 
of stocks with similar correlations, (corresponding to Market Sectors). 

Using the correlation matrix from the "one-factor model" and the "market 
plus sectors model" , we can construct correlated time series using the Cholesky 
decomposition A of a correlation matrix C = AA T . We can then generate finite 
correlated time series of length T, 

x it = AijVjt t = l,...,T (6) 
j 



where yjt is a random Gaussian variable with mean zero and variance 1 at 
time t. Using Eqn. [2] we can then construct a correlation matrix using the 
simulated time series. The finite size of the time series introduces 'noise' into 
the system and hence empirical correlations will vary from sample to sample. 
This 'noise' could be reduced through the use of longer simulated time series 
or through averaging over a large number of time series. 

In order to compare the eigenvectors from each of the Model Correlation 
matrices to that constructed from the equity returns time series, we use the 
Inverse Participation Ratio (IPR) |4"f29~] . The IPR allows quantification of the 
number of components that participate significantly in each eigenvector and 
tells us more about the level and nature of deviation from RMT. The IPR of 
the eigenvector u k is given by I k = J2iLi ( u i) an d allows us to compute the 
inverse of the number of eigenvector components that contribute significantly 
to each eigenvector. 



3 Data 



In order to study the dynamics of the empirical correlation matrix over time, 
we analyse two different data sets. The first data set comprises the 384 equi- 
ties of the Standard & Poors (S&P) 500 where full price data is available from 
January 1996 to August 2007 resulting in 2938 daily returns. The S&P 500 
is an index consisting of 500 large capitalisation equities, which are predomi- 
nantly from the US. In order to demonstrate that our results are not market 
specific, however, we examine a second data set, made up of the 49 equities of 
the Dow Jones Euro Stoxx 50 where full price data is available from January 
2001 to August 2007 resulting in 1619 daily returns. The Dow Jones Euro 
Stoxx 50 is a stock index of Eurozone equities designed to provide a blue-chip 
representation of supersector leaders in the Eurozone. 
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4 Results 



We analyse the eigenvalue dynamics of the correlation matrix of a subset of 
100 S&P equities, chosen randomly, using a sliding window of 200 days. This 
subsector was chosen such that Q = = 2, thus ensuring that the data 
would be close to non-stationary in each sliding window. Figure Ufa) shows 
broadly similar sample dynamics from the 5 , 15 th and 25 th largest eigenvalues 
over each of these sliding windows. The sum of the 80 smallest eigenvalues 
are shown in Figure d](b), while the dynamics of the largest eigenvalue is 
displayed in Figure [0(c). The repulsion between the largest eigenvalue and 
the small eigenvalues is evident here, (comparing [H(b) and [U^c)), with the 
dynamics of the small eigenvalues contrary to those of the largest eigenvalue. 
As noted earlier (Section 12. ip . this is a consequence of the fact that the trace 
of the correlation matrix must remain constant under transformations and any 
change in the largest eigenvalue must be reflected by a change in one or more 
of the other eigenvalues. Similar results were obtained for different subsets of 
the S&P and also for the members of the Dow Jones Euro Stoxx 50. 




1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 




1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 

Time (Years) 



Fig. 1. Time Evolution of (a) Three small eigenvalues (b) Sum of the 80 smallest 
eigenvalues (c) The largest eigenvalue 
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4-1 Normalised Eigenvalue Dynamics 

Using normalised eigenvalues as described above, (Eqn. EJ), we performed a 
number of experiments to investigate the dynamics of a set of small eigenvalues 
versus the largest eigenvalue. The various experiments are described below: 

(1) As in Section HI the dynamics for the same subset of 100 equities are 
analysed using a sliding window of 200 days. The normalisation is carried 
out using the mean and standard deviation of each of the eigenvalues over 
the entire time-period. Figure [2(a) shows the value of the S&P index from 
1997 to mid-2007. 

The normalised largest eigenvalue is shown in Figure [21(b) together 
with the average of the 80 normalised small eigenvalues. The compen- 
satory dynamics mentioned earlier are shown more clearly here, with the 
largest and average of the smallest 80 eigenvalues having opposite move- 
ments. The normalised eigenvalues for the entire eigenvalue spectrum are 
shown in Figure 12(c), where the colour indicates the number of standard 
deviations from the time average for each of the eigenvalues over time. As 
shown, there is very little to differentiate the dynamics of the 80 — 90 or so 
smallest eigenvalues. In contrast, the behaviour of the largest eigenvalue 
is clearly opposite to that of the smaller eigenvalues. However, from the 
90 th and subsequent eigenvalue there is a marked change in the behaviour, 
(Figure [2(d)), and the eigenvalue dynamics are distinctly different. This 
may correspond to the area outside the "Random Bulk" in RMT. Simi- 
lar to [24.25J , we also find evidence of an increase/decrease in the largest 
eigenvalue with respect to 'drawdowns'/'draw-ups'. Additionally, we find 
the highlighted compensatory dynamics of the small eigenvalues. These 
results were tested for various time windows and normalisation periods, 
with smaller windows and normalisation periods found to capture and 
emphasise additional features. 

(2) To demonstrate the above result for a different level of granularity, we 
chose 50 equities randomly with a time window of 500 days, giving Q = 

= 10. The results obtained, (Figure EJ), are in keeping with those for 
Q = 2 earlier, with a broad-band increase (decrease) of the 40 smallest 
eigenvalues concurrent to a decrease (increase) of the largest eigenvalue. 

(3) The previous examples used random subsets of the S&P universe in order 
to keep Q = ^ as large as possible. To demonstrate that the above results 
were not sampling artifacts, we also looked at the full sample of 384 
equities, (that survived the entire 11 year period), with a time window 
of 500 days (Q = 1.30). The results, as shown in Figured! are similar to 
those above, with the majority of the small eigenvalues compensating for 
changes in the large eigenvalue. As indicated previously, however, there 
is a small band of large eigenvalues, for which behaviour is different to 
that of both the band of small eigenvalues and the largest eigenvalue. 
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Fig. 2. (a) S&P Index (b) Normalised Largest Eigenvalue vs. Average of 80 smallest 
normalised eigenvalues (c) All Normalised Eigenvalues (d) Largest 12 Normalised 
Eigenvalues 
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Fig. 3. (a) Normalised Largest Eigenvalue vs. Average of 40 smallest normalised 
eigenvalues (b) All Normalised Eigenvalues 
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Fig. 4. (a) Normalised Largest Eigenvalue vs. Average of 325 smallest normalised 
eigenvalues (b) All Normalised Eigenvalues 

(4) All examples discussed so far have focused on the universe of equities 
from the S&P 500 that have survived since 1997. To ensure that the re- 
sults obtained were not exclusive to the S&P 500, we also applied the 
same technique to the 49 equities of the EuroStoxx 50 index that sur- 
vived from January 2001 to August 2007. The sliding window used was 
200 days, such that Q = 4.082. The results, (Figure E]), were again similar 
to before, with a wide band of small eigenvalues "responding to" move- 
ments in the largest eigenvalues. In this case, the band of deviating large 
eigenvalues (ie. those which correspond to the area outside the "Random 
Bulk" in RMT), (Figure E^d)), is not as marked as in the previous exam- 
ple. This effectively implies that equities in this index are dominated by 
"the Market" . 



4-2 Model Correlation Matrix 



The results, described, demonstrate that the time dependent dynamics of the 
small eigenvalues of the correlation matrix of stock returns move counter to 
those of the largest eigenvalue. Here, we show how a simple one-factor model, 
Section 12.21 of the correlation structure reproduces much of this behaviour. 
Furthermore, we show how additional features can be captured by including 
perturbations in this model, essentially a "market plus sectors" model, Sec- 
tion [231 [2211291130] . 

In order to compare the empirical results, Section HI to those of the single fac- 
tor model, we first constructed a correlation matrix where each non-diagonal 
element was equal to the average correlation of the empirical matrix in each 
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Fig. 5. (a) EuroStoxx 50 Index, Jan 2001 - Aug 2007 (b) Normalised Largest Eigen- 
value vs. Average of 40 smallest normalised eigenvalues for EuroStoxx 50 (c) All 
Normalised Eigenvalues (d) The 9 Largest Normalised Eigenvalues 

sliding window. We then calculated the eigenvalues of this matrix over each 
sliding window and normalised these as before, (Section [2} . The dynamics of 
the largest eigenvalue for the single-factor model are displayed in Figure Efa) 
for the EuroStoxx 50 index using a sliding window of 200 days. As can be seen, 
the main features of the dynamics are in agreement with those of Figure [5] for 
the empirical data. The dynamics of the remaining eigenvalues, shown in fig- 
ure EIc), are found to be within a narrow range, with any change in time due 
to compensation for change in the largest eigenvalue. 

For the 'market plus sectors' model, we introducted perturbations with two 
groups of 5 stocks having correlation p — 0.15 and p Q + 0.15, with the average 
correlation at each time window remaining the same. The largest eigenvalue 
dynamics for the 'market plus sectors' model is shown m^b), with the dy- 
namics almost identical to those for the largest eigenvalue of the 'one-factor' 
model, (any differences are due to random fluctuations induced in the Cholesky 
decomposition, Section [2.3p . However, looking at the remaining eigenvalues, 
E{cl), we see that a number are found to be deviate significantly from the bulk. 
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These deviations are due to the additional sectorial information included in 
the market plus sectors model. This agrees with previous results, [HE], where 
small eigenvalues were found to deviate from the random bulk, in addition to 
large deviating eigenvalues containing sectorial information. 




(b) i6 
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2002 2003 2004 2005 2006 2007 

Time 



Fig. 6. (a) Largest Eigenavalue, one factor model (b) Largest Eigenavalue, Market 
plus sectors model (c) Eigenvalues 1-48, one factor model (d) Eigenvalues 1-48, 
Market plus sectors model 

To examine properties of the eigenvector components themselves, we use the 
Inverse Participation Ratio (IPR). Figure [7(a), displays the IPR found using 
the empirical data from the Eurostoxx 50. This has similar characteristics to 
those found for different indices, [I], with the IPR for the largest eigenvalue 
much smaller than the mean, a large IPR corresponding to sectorial informa- 
tion in the 2 nd or 3 rd largest eigenvalues and the IPR for the small eigenvalues 
raised. For the single factor model, we created a synthetic correlation matrix 
using Eqn. with average correlation (0.204) equal to that of the Euro 
Stoxx 50 over the time period studied. As shown in Figure UJb), the IPR 
retains some of the features found for empirical data, [4,29j, with that cor- 
responding to the largest eigenvector having a much smaller value than the 
mean. This corresponds to an eigenvector to which many stocks contribute, 
(effectively the market eigenvector), [4"f29] . 

In an attempt to include additional empirical features, such as the band of 
deviating large eigenvalues between the bulk and the largest eigenvalue, we also 
considered a perturbation, with two groups of stocks having correlation p — 
0.15 and p + 0.15, and the same average correlation. In this case, Figure UJb), 
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additional features are found, with larger IPR for both smallest and second 
largest eigenvalue. This agrees with |4] where, for empirical data, the group 
structure resulted in a number of small and large eigenvalues with larger IPR 
than that of the bulk of eigenvalues. These large eigenvalues were shown, |4"|IT0"] , 
to be associated with correlation information related to the group structure. 



0.5 




O 5 10 15 20 25 30 35 40 45 50 




02 I 1 1 1 1 1 1 1 1 1 1 

5 10 15 20 25 30 35 40 45 50 



Fig. 7. (a) IPR Empirical Data (b) IPR Simulated Data 



4-3 Drawdown Analysis 



As described, [21], and discussed above, drawdowns (periods of large negative 
returns) tend to reflect an increase of one eigenstate of the cross-correlation 
matrix, with the opposite occurring during drawups. In this section, we at- 
tempt to characterise the market according to the relative size of both the 
largest and small eigenvalues using stocks of the Dow Jones Euro Stoxx 50. 

Returns, correlation matrix and eigenvalue spectrum time series for overlap- 
ping windows of 200 days were calculated and normalised using the mean and 
standard deviation over the entire series, (Eqn. [5]). Representing normalised 
eigenvalues in terms of standard deviation units (SDU) allows partitioning 
according to the magnitude of the eigenvalues. Table [1] shows the average re- 
turn of the index, (during periods when the largest eigenvalue is ±1 SDU), 
for both the largest eigenvalue and the average of the normalised 40 smallest 
eigenvalues. 

The results, Table [TJ demonstrate that when the largest eigenvalue is > 1 SDU, 
the average index return over a 200 day period is found to be —16.80%. When 
it is small (< —1 SDU), the average index return is 18.46%. Hence, the largest 
eigenvalue can be used to characterise markets, with the eigenvalue peaking 
during periods of negative returns (Drawdowns) and bottoming out when the 
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market is rising (Drawup). For the average of the 40 smallest eigenvalues, 
the partition also reflected drawdowns and drawups, but with opposite signs. 
This indicates that information about the correlation dynamics of financial 
time series is visible in both the lower and upper eigenstates, in agreement 
with the results found by [T3|26] for both synthetic data and EEG seizure 
data. 

Eigenvalues No. Std Index Return 



Large 
Average 40 Smallest 



<-l 18.46% 
>1 -16.80% 



<-l -23.yu7o 
>1 19.32% 



Table 1 

Drawdown/Drawup analysis, average index returns for eigenvalue partitions in SDU. 



5 Conclusions 



The correlation structure of multivariate financial time series was studied by 
investigation of the eigenvalue spectrum of the equal-time cross-correlation 
matrix. By filtering the correlation matrix through the use of a sliding window 
we have examined behaviour of the largest eigenvalue over time. As shown, 
Figures (T5]-[5]), the largest eigenvalue moves counter to that of a band of small 
eigenvalues, due to eigenvalue repulsion. A decrease in the largest eigenvalue, 
with a corresponding increase in the small eigenvalues, corresponds to a re- 
distribution of the correlation structure across more dimensions of the vector 
space spanned by the correlation matrix. Hence, additional eigenvalues are 
needed to explain the correlation structure in the data. Conversely, when the 
correlation structure is dominated by a smaller number of factors (eg. the 
"single-factor model" of equity returns), the number of eigenvalues needed to 
describe the correlation structure in the data is reduced. In the context of pre- 
vious work, P4"f25] . this means that fewer eigenvalues are needed to describe 
the correlation structure of 'drawdowns' than that of 'draw-ups'. 

By introducing a simple 'one-factor model' of the correlation in the system 
(Section 14. 2j) . we were able to reproduce the main results of the empirical 
study. The compensatory dynamics, described, were clearly seen for a correla- 
tion matrix with all elements equal to the average of the empirical correlation 
matrix. The model was then adapted, by the addition of pertubations to the 
correlations, with the average correlation remaining unchanged. This 'mar- 
kets plus sectors' type model was able to reproduce additional features of 
the empirical correlation matrix with eigenvalues deviating from below and 
above the bulk. The Inverse Participation Ratio of the "markets plus sectors" 
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model was also shown to have group characteristics typically associated with 
known Industrial Sectors, with a larger than average value for the smallest 
eigenvalue and for the second largest eigenvalue. Through a partition of the 
time-normalised eigenvalues, it was then shown quantitatively that the largest 
eigenvalue is greatest /smallest during drawdowns/drawups, and vice versa for 
the small eigenvalues. 

Suggested future work includes a more detailed study of the relationship be- 
tween the direction of the market and magnitude of the eigenvalues of the 
correlation matrix. Studying the multiscaled correlation dynamics over dif- 
ferent granularities may shed some light on the different collective behaviour 
of traders with different strategies and time horizons. Additional analysis of 
high frequency data may also be useful in the characterisation of correlation 
dynamics, especially prior to market crashes. Study of the possible relation- 
ship between the dynamics of the small eigenvalues and additional correlation 
information which, according to some authors, |19f2Uf21f22f23f26f27j . may be 
hidden in the part of the eigenvalue spectrum normally classifed as noise could 
be achieved through analysis of the relative dynamics of the small and large 
eigenvalues at times of extreme volatility, (such as during market crashes). 
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